The Shape of a Ponytail and the Statistical Physics of Hair Fiber Bundles 
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A general continuum theory for the distribution of hairs in a bundle is developed, treating indi- 
vidual fibers as elastic filaments with random intrinsic curvatures. Applying this formalism to the 
iconic problem of the ponytail, the combined effects of bending elasticity, gravity, and orient at ional 
disorder are recast as a differential equation for the envelope of the bundle, in which the compress- 
ibility enters through an 'equation of state'. From this, we identify the balance of forces in various 
regions of the ponytail, extract a remarkably simple equation of state from laboratory measurements 
of human ponytails, and relate the pressure to the measured random curvatures of individual hairs. 



PACS numbers: 87.19.R-, 05.45.-a, 46.65. +g 

One of the most familiar features of a bundle of hair 
such as a ponytail is its 'body' or 'volume.' Close ex- 
amination reveals that this property arises in a sub- 
tle way from the stiffness and shapes of the individual 
fibers, whose meandering paths through the bundle pro- 
duce many collisions with other hairs (Fig. [T]). These 
meanderings are in part a consequence of the contacts 
themselves, but hairs also have an intrinsic waviness or 
curliness [U [2] . Such curvatures may be generated dur- 
ing growth, and vary with ethnicity. They are clearly also 
modified by chemical, thermal, and mechanical forces, as 
in the 'water wave' treatment, or a 'perm' [1 . 

From Leonardo to the Brothers Grimm our fascination 
with hair has endured in art and science [3, 4 . Yet, 
we still do not have an answer to perhaps the simplest 
question that captures the competing effects of filament 
elasticity, gravity, and mutual interactions: What is the 
shape of a ponytail? Note that the average human has 
^ 10^ head hairs, so if even a modest fraction is gathered 
into a ponytail, the number involved is enormous: this is 
a problem in statistical physics. 

Here we propose a theory for the ponytail shape on the 
basis of a continuum theory for the spatial distribution of 
hairs in a bundle. Their random curvatures give rise to a 
swelling pressure characterized by an 'equation of state' 
(EOS) of hair, a concept first introduced semi-empirically 
by van Wyk in 1946 in relation to the compressibility of 
wool O [6], and explored recently for two-dimensional 
randomly- curved fibers by Beckrich et al. [7 . We ex- 
ploit the predominantly vertical alignment of hairs and 
axisymmetry to justify a number of approximations that 
render the problem analytically tractable, and thereby 
reduce the many-body problem to a one-body problem 
for the ponytail envelope. This shows how the EOS mod- 
ifies the envelope shape from that of a single hair bent 
by gravity, a classic problem in elasticity [2 . 

In parallel with the theoretical development, we mea- 
sured the shapes of ponytails made from commercially 



available hair 'switches' [8], and of their component 
fibers. Typical human hair has an elliptical cross section 
and a distribution of major axis diameters 40 < (i < 140 
Ilia. We found d = 79 ± 16 /im for a random sam- 
ple from the switches. Hair has an average density of 

1.3g/cm^ [1 , and a linear mass density A ^ 65 jig/ cm 
(in more amusing units, 6.5 g/km). Though its inter- 
nal microstructure is complex, the bend and twist mod- 
uli of hairs [l] are consistent with those of a homoge- 
neous incompressible material with a nylon-like modu- 
lus E ^ A GPa. On the centimeter scale classical fil- 
ament elasticity holds, with a bending modulus A = 
Endt^/M ^ 8 X 10"^ Nm^ The quantities A and A and 
the acceleration of gravity g combine to form the length 
I = {A/Xg)^^^ ^ 5 cm on which gravity bends a hair [2]. 

Individual hairs display a range of shapes (Fig. ^) 
which we have quantified by high-resolution stereo- 
scopic imaging [9 . Both the mean squared curvature 




FIG. 1: (color online) A ponytail. (a) Rot at ionally- averaged 
image of a switch oi N ^ 9500 fibers, approximately 25 cm 
long. Coordinate system for envelope shape R{s) in terms of 
arc length s{z). (b) Meanderings of hairs near ponytail edge. 
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FIG. 2: (color online) Geometry of hairs, (a) Representative projections of hair contours, thickened for clarity, (b) Mean squared 
radial excursion and curvature as functions of arc length, from processing stereoscopic image pairs [9 . The reconstructed arc 
lengths cluster tightly around 24.50 it 0.05cm (inset histogram), demonstrating the accuracy of the image processing and 
analysis methodology. Error bars are standard errors from ensemble averaging (N =115 fibers in total). 



{kappao{s)'^) and the radial excursion (ro(s)^) increase 
with arclength s measured from the top of the switch 
(Fig. Whilst some of this is undoubtedly due to 

gravity (recall £ ^ 5 cm), the major part is intrinsic, as 
we have verified by examining inverted hairs. This is in 
part due to the preparation process: after washing and 
rinsing, the hairs in a drying ponytail pass through a glass 
transition with decreasing humidity [1 , locking in the in- 
trinsic curvature [lOj, which is naturally reduced in the 
vicinity of the clamp due to confinement by neighboring 
filaments. Although this is something of a complication 
when it comes to interpreting the results, we must re- 
gard it as an essential feature of hair switches and pony- 
tails comprised of real fibers. For later reference, the 
length- wise averages are (/^q) = 0.15 ± 0.01 cm~^ and 
(rg) = 0.80 ± 0.05 cm^. 

Figure [S^i shows measured profiles of radius R{z) vs 
distance z below the clamp for four separate switches of 
length L ^ 25 cm. Each profile has been obtained from 
the rotational average of five images as in Fig. [l^ viewed 
from angles 72° apart. The switch profile shows quite 
good reproducibility and is well modeled and explained 
by the theory we now describe. Our starting point for 
the continuum theory is to introduce the fiber length 
density p{r) (the number of fibers per unit area inter- 
secting a plane perpendicular to the fibers) and the mean 
fiber tangent vector t(r), the local average of unit vec- 
tors along the fibers. The latter is a meaningful quantity 
when the fiber orientation remains coherent over length 
scales much larger than the mean fiber spacing p~^l^ . 
Here the fibers are indeed well- aligned, with t = |t| ~ 1, 
unlike in non- woven fabrics [11]. In the absence of fiber 
ends in the bulk these continuum fields obey a continu- 
ity equation V • (pt) =0. The analogy to the continu- 
ity equation of fluid mechanics mathematizes the remark 
made by Leonardo at the beginning of the 16th century, 
that hair resembles fluid streamlines [3], an observation 



which has been exploited in more recent times to aid com- 
puter animation [4 . For later use we also define the local 
packing fraction ^ = npdP/A. We propose the energy of 
an axisymmetric fiber bundle is 



£[p,t] = jc 



(fvp -Ak^ + ip{r) + (u) 



(1) 



where hi = |(t • V)t| is the curvature field. The terms in 
([T]) are the elastic energy of mean curvature, the external 
{e.g. gravitational) potential (f^ and a fiber confinement 
energy per unit length (u) that aggregates all terms in- 
volving disorder, such as contacts and natural curvatures. 
Without axisymmetry, ([T]) should include terms arising 
from the torsion of t. As in density functional theory 
^12] , we suppose that (u) is some local function of p. 
Minimization of ([T]) provides a variational principle for 
the bundle shape and the distribution of fibers. When 
recast as mechanical force balance we make contact with 
the EOS, and identify P{p) = p^d{u) /dp as the pressure. 

To address the specific problem of ponytail shapes we 
now introduce models which allow for largely analyti- 
cal calculations. With axisymmetry, an integrated form 
of the continuity equation is 27rr/)sin^ = —dn/dz and 
27TrpcosO = dn/dr^ where n(r, z) is the number of fibers 
within a radius r at depth z and 6 is the angle the tangent 
vector makes to the vertical. (n(r, z) plays the role of the 
stream function in fluid mechanics.) We insert this into 
Eq. ([T]), and use a trial uniform radial density function 
with s elf- similar form ^ n{r,z) = N[r / R{z)]'^ where N is 
the total number of fibers and R{z) is the ponytail radius 
(Fig. [l^). In practice it is more convenient to use R{s), 
where s{z) is the arclength from the clamp. If L is the 
fully-extended hair length and (p = Xgz the gravitational 
potential energy, then to second order in Rg = dR/ds, 
neglecting a small splay term, one finds 
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FIG. 3: (color online) Analysis of ponytail shapes, (a) Measured hair switch profiles (thin black lines), compared to the 
prediction of Eq. (jsj with n(i?) as given (solid blue line), with n(i?) = (dashed blue line), and with the van Wyk EOS 
(dashed red line) [51113] (b) Dimensionless swelling pressure Il{R) from cubic and quartic fits to the measured profiles, using 
procedure outlined in the text, (c) The magnitude of the four terms in Eq. Jsl for the calculated profile (solid blue line) in (a). 



where (u) depends on p = N/{t:R^). The problem is 
now mapped to an equivalent single fiber hanging under 
gravity in the presence of a radial force field derived from 
{u). The uniform distribution in the trial density func- 
tion yields renormalized material properties A = A/2 
and A = A/2. Minimizing Eq. ([2| leads to 

l^Rssss -{L- s)Rss ^Rs- n(i?) = (3) 

where U{R) = 4£^P/ApR = -(2^ /A) d{u) / dR. We 
term this the ponytail shape equation. It describes a force 
balance on a length element of the notionally equivalent 
single fiber as the sum of four dimensionless terms which 
are, respectively, an elastic restoring force, a 'string ten- 
sion' contribution, a weight term, and a radial swelling 
force corresponding to a pressure gradient P/ R per unit 
fiber density. The ratio Ra = L/l we shall term the Ra- 
punzel number^ since it is a dimensionless measure of the 
ponytail length. When the ponytail hangs from a cir- 
cular clamp of radius Rc^ the boundary conditions are 
i?(0) = Rc and i?s(0) = tan^c where 9c is the 'launch' 
angle of the outermost fibers emerging from the clamp. 
At the free bottom of the ponytail the boundary condi- 
tions are Rss{L) = Rsss{L) = 0. To the order at which 
we are working, (|3| is supplemented hy Zs — I — Rl/2 to 
give the parametric ponytail shape {z{s), R{s)). 

Fitting the above theory to the experimental ponytail 
profiles in Figure [3^ reveals a remarkably simple form 
for the pressure U{R). While the Ml Eq. ^ can m 
principle be used to determine Il{R) from the profiles, 
the extraction of high-order derivatives from such data 
is notoriously problematic. We notice though that, away 
from the clamp, Rssss is likely to be subdominant to 
the other terms in Eq. Q and therefore we can neglect 
this elastic term and approximate Il{R) — Rg — {L — 
s)Rss^ where the right-hand-side is obtained by a low- 
order polynomial fit to the data. Figure ^jp shows that 



in this region the EOS is accurately represented by 

U{R) = Uo{l-R/Rn , (4) 

with Ho ~ 0.85 and R* ^ 6 cm. Inserting this into Eq. (|3| 
and now including the elastic term recovers the solid blue 
line in Fig. [3|i, in excellent agreement with the data (by 
contrast the van Wyk EOS simply cannot be made to 
fit the data [l3]). In making these calculations we use 
Rc ~ 4 mm and Oc ^ 17°, obtained from measurements 
near the clamp. The starting radius Rc corresponds to 
(j) ^ 0.95, consistent with the near close packing of the 
fibers, whilst the starting angle 6c is presumably governed 
by the method of clamping (in our case a rubber band 
wrapped several times around the top of the switch). 

Figure|3]3 shows the magnitudes of the terms in Eq. 
In the region near the clamp {s < 2 cm), elasticity and 
pressure balance, but for the most part the dominant 
balance is between weight and pressure, justifying our 
claim that the elastic term is subdominant away from 
the clamp. The blue dashed line in Fig. |3^ is the profile 
for Il{R) = 0. Since Ra ~ 5 is quite large, this shape is 
dominated by gravity. Comparing the dashed and solid 
blue lines in Fig. |3^ highlights again the dominant role 
played by the swelling pressure in determining the shape. 

Given n(i?), the shape of any ponytail can be pre- 
dicted. Thus we are led to a kind of experimentum cru- 
cis^ shown in Fig. [4j in which the predictions of Eq. (|3| 
are compared to the profiles of progressively cut hair 
switches. The agreement is very good. We observe em- 
pirically that the launch angle Oc is remarkably constant 
at 17°, only decreasing to 16° for the shortest hair switch 
(all calculations used Oc = 17°). The calculated profiles 
show a modest compaction on increasing length, while 
the experimental profiles almost completely collapse on 
top of one another. This is not an effect of plasticity [11] 
since the switches are compressed in the cutting stage. 
The predicted profiles can similarly collapse (Fig.jlj^), by 
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FIG. 4: Trimming a ponytail. (a) Superimposed rotationally- 
averaged outlines of four hair switches, cut down from 25 cm 
in steps of 5 cm. (b) Predicted profiles from Eq. (Is]) with 
U{R) = no(l - (c) Predicted profiles with n(S, s) = 

no(l - R/R*) X 2s/ L* where L* = 25cm. 



inter-fiber attractions, such as the capillary forces acting 
on wet hair [TF or a paintbrush. 

The program laid out here extends some central 
paradigms in statistical physics to the enchanting prob- 
lem of ponytail shapes. The remarkably simple equation 
of state we have found, along with the systematic vari- 
ation of intrinsic curvature along fibers, may open the 
way to understanding a wide range of hair and fur ge- 
ometries. It is also of interest to extend the analysis to 
the dynamics of fiber bundles, epitomized by the 'swing' 
of a ponytail [16 , where the notion of an equivalent single 
fiber may again prove fruitful. 

We thank A. Avery, M.E. Gates, and A.S. Ferrante for 
helpful discussions. This work was supported in part by 
the Schlumberger Ghair Fund. 



allowing for an additional length dependence reflecting 
the gradient in the intrinsic fiber properties (Fig. [2]3). 

How are we to interpret the EOS recovered by this 
analysis? We propose that it can be associated with the 
intrinsic curvatures of the filaments. Let us imagine that 
the effects of collisions with neighboring fibers can be 
captured by a tube model. Specifically, consider a helical 
fiber \]A of radius ao confined within a cylinder of ra- 
dius a < ao, for which (u) ^ A(/^q)(1 — a/ao)^/2 (this is 
not exact but is quite accurate). Matching this to (u) = 
{A/2i^) n(i?) dR^ obtained by integrating the expres- 
sion below Eq. ([3|, and inserting our empirical Il{R) gives 
a/ao ^ l-a^aR/R"" where a = ^/[l^WJW{K^ ^ OA. 
If we additionally suppose ^ (r^) then for instance at 
the clamp {R = Re) the confining tube radius a ^ 6 mm. 
This seems to be in reasonable accord with observations, 
(see e. g. Fig. ^p)- In the tube model, the pressure thus 
arises from increasing confinement of the fibers but they 
are still far from being completely straightened, even at 
high compression. This is perhaps not surprising consid- 
ering the role that must eventually be played by friction. 

Of the existing EOS theories, that of van Wyk does not 
fit our data (red dashed line, Fig.[3|i), nor does it have an 
explicit link to the random curvatures. The fiber-collision 
model [7] can be extended to three dimensions, but the 
link to the underlying statistical properties becomes very 
unwieldy. More importantly, in that model fiber excur- 
sions are limited by nearest-neighbor collisions. This is 
not necessarily the case in three dimensions, and in fact 
is not supported by our data. Hence the microscopic link 
between fiber confinement and packing fraction remains 
an important open problem. Interestingly both the tube 
model and the collision model predict that the pressure 
remains finite on approach to the close packing limit, in 
marked contrast to thermal systems of hard particles. 
Thus, a bundle can be collapsed by sufficiently strong 
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